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ABSTRACT 



^—i , Aims. The a- and 7-effects, which are responsible for the generation and turbulent pumping of large scale magnetic fields, 
fNj ■ respectively, due to passive advection by convection are determined in the rapid rotation regime corresponding to the deep 
layers of the solar convection zone. 
I Methods. A 3D rectangular local model is used for solving the full set of MIfD equations in order to compute the electromotive 
' force (emf), £ — u x b, generated by the interaction of imposed weak gradient-free magnetic fields and turbulent convection 
(-H _ with varying rotational infiuence and latitude. By expanding the emf in terms of the mean magnetic field. Si = aijBj, all 
^ nine components of aij are computed. The diagonal elements of Uij describe the a-effect, whereas the off-diagonals represent 
magnetic pumping. The latter is essentially the advection of magnetic fields by means other than the underlying large-scale 
velocity field. Comparisons are made to analytical expressions of the coefficients derived under the first-order smoothing 
approximation (FOSA). 

] Results. In the rapid rotation regime the latitudinal dependence of the a-components responsible for the generation of the 
• • ■ azimuthal and radial fields does not exhibit a peak at the poles, as is the case for slow rotation, but at a latitude of about 
30° . The magnetic pumping is predominantly radially down- and latitudinally equatorward as in earlier studies. The numerical 
results compare surprisingly well with analytical expressions derived under first-order smoothing, although the present 
^ calculations are expected to lie near the limits of the validity range of FOSA. 

_ _ I 
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1. Introduction 

According to dynamo theory small-scale helical turbulence 
caused by the convective instability in combination with 
global rotation, uniform or differential, is the source of the 
large-scale magnetic structures seen in the Sun and in var- 
ious other late-type stars (e.g. Parker 0*9551 MofFatt [19781 
Krause & Radler ITg8D|) . The small scales enter the evolu- 
tion equation of the mean magnetic field via the turbulent 
electromotive force £ = u x b. li the mean magnetic field 
B varies slowly in time and space, the emf can be repre- 
sented in terms of B and its gradients 

0~B ■ 

Si = aijBj + bi-jk^r^ + ... , (1) 
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where aij and bijk are in the general case tensors contain- 
ing the transport coefficients, and the dots indicate that 
higher order derivatives can be taken into account. The 
challenge is to derive the tensors Uij and bijk which, in 
general, cannot be done from first principles due to the 
lack of a comprehensive theory of convective turbulence. 

In the kinematic regime where the magnetic energy 
is negligible in comparison to the kinetic energy the 
most simple approximation is to neglect all correlations 
higher than second order in the fluctuations. This is of- 
ten called the first order smoothing (FOSA), quasilinear, 
or the second order correlation approximation (SOCA). 
Considering the simple case of isotropic turbulence in the 
high-conductivity limit as an example where the transport 
coeffients can be computed analytically, the tensor aij re- 
duces into a single scalar (Steenbeck & Krause [l969|l 

a^-^TcuT^, (2) 
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where is the correlation time of the turbulence, uj = 
V X It the vorticity, and w • u the kinetic helicity. We 
note that in the anisotropic case the trace of the tensor 
Oij, instead of any individual component of the a-effect, is 
proportional to the kinetic helicity (e.g. Radler fl980l see 
also Sect. EH). 

Despite the limited applicability of Eq. numeri- 
cal models of magnetoconvection have shown that there is 
at least a qualitative agreement between the a-effect and 
the negative of the kinetic helicity also in the anisotropic 
case (Brandenburg et al. 119901 Ossendrijver et al. IBUOTI 
hereafter Paper I). Furthermore, using local time-distance 
helioseismology, Gizon & Duvall ( 2003 ) studied the corre- 
lation of the curl and divergence of the horizontal compo- 
nents of solar surface flows as a proxy for the helicity (see 
also Riidiger et al. .1999^ and found that the correlation 
follows a cos^ latitude profile for the observed latitude 
range. Using local numerical convection models Egorov et 
al. ((^004) were able to reproduce the cos 6* latitude profile 
and further confirmed the validity of the correlation as a 
valid tracer of the helicity. In a recent study, Kapyla et al. 
H2004|) found a similar latitude dependence for moderate 
rotation, i.e. up to Coriolis number Co — 2fll/u < 4, in 
agreement with the results of Egorov et al. For more rapid 
rotation, Co ~ 10, however, the same authors found that 
the latitude distribution of the volume average over the 
convectively unstable region resembles more a sin 6 profile 
for latitudes higher than 15°. As our working hypothesis, 
we consider the relation (0) to trace the a-effect, which 
raises the question of the latitude dependence of this ef- 
fect in the rapid rotation regime. 

In the present study we examine the a-effect by means 
of local numerical modelling of small volumes of a star 
at different latitudes. In distinction to other studies (e.g. 
Brandenburg et al. [lp90, Ossendrijver et al. IMTTl I^HH^ 
hereafter Paper II), we extend the parameter range in the 
rapid rotation regime to values typical for the base of the 
solar convection zone, where the Coriolis number is of the 
order of 10 or larger (see e.g. Kiiker et al. 119931 Kapyla et 
a,].l2005h|). 

As a secondary objective we make an attempt to test 
the validity of FOSA. A sufficient, but not a necessary, 
requirement for the applicability of FOSA is that either 
the Reynolds or the Strouhal number is small 



iml — -jUfcfTc 1 



min(Rm, St) <^ 1 , 



(3) 



where u and I are the typical velocity and length scales, 
fcf the wavenumber of the energy carrying scale which 
is related with the correlation lenght via Ic = 27r/fcf. In 
the high-conductivity limit relevant for Eq. (O, where 
rjTc/ll <C 1, the requirement Q reduces simply to St ^ 
1 which, however, is at best marginally satisfied (e.g. 
Brandenburg et al. 120041 Brandenburg & Subramanian 
I2005bl Kapyla et aL lUnSj) . 

We derive analytical formulae for all of the ay- 
coefficients in the high-conductivity limit. For the different 
components of the a-effect we find expressions analogous 



to Eq. ^ in the anisotropic case. We then compare the 
analytical result and its numerical equivalent using the 
correlation time as a free parameter. The correlation time 
allows us, at least in principle, to compute the Strouhal 
number. 

The remainder of the paper is organised as follows: in 
Sect. 12 the numerical convection model is described. In 
Sect. O the problem and research methods are stated. In 
Sects. 0] and |31 we give the results and conclusions, respec- 
tively. 



2. The convection model 

2.1. Basic equations 

The convection model is the same as that presented in 
Kapyla et al. ( 2004 , 2005b) . The computational domain is 
a rectangular box at latitude O in the southern hemisphere 
of a star. The coordinate system of the model is chosen 
such that X, y, and z correspond to (—0, —r) in spherical 
coordinates, where 9 = 90° — 6 is the colatitude. The 
angular velocity as a function of latitude is given by = 
r2(cos9ea, — sinOe^). The box has horizontal dimensions 
= Ly = 4, and = 2 in the vertical direction in 
units of the depth of the convectively unstable layer d. 
We set (zq, zi, Z2, 2:3) = (—0.15, 0, 1, 1.85), where zq and Z3 
correspond to the top and bottom boundaries of the box, 
respectively, whereas zi and Z2 give the positions of the 
upper and lower boundaries of the convectively unstable 
layer. We solve a system of MHD-equations 



— — = u X B — tii.j,qJ 



dA 

'dt 
91np 

dt 
du 

'dt 



de 
dt 



-(m • V) ln/7 + V • M , 



(4) 
(5) 



1. 



-(u ■ V)u Vp -2nxu + g + 

P 

+ -J xB + -\7 ■ {2vpS) , (6) 
P P 

-(«• V)e-^(V-«) + ^j2 + 
P P 

- rcooi , (7) 



where A is the magnetic vector potential, u the velocity, 
B = V X A + Bq the magnetic field which is the sum 
of the fiuctuating and (constant) imposed contributions, 
J = V X B / fj,Q the current density, p the mass density, p 
the pressure, g = gSz the constant gravity, the unit vec- 



tor in the vertical direction. Si 



2 



the strain tensor, and e — cyT the internal energy per 
unit mass, po is the vacuum permeability, 77 and v the con- 
stant magnetic diffusivity and kinematic viscosity. Fcond 
describes thermal conduction (see the next subsection). 
The term Fvisc describes viscous heating via 



2i/Si 



duj 
dxi 



(8) 
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The uppermost layer of the box is cooled according to 



1 



cool 



^cool 



fiz){e - eo) , 



(9) 



where tcooi is a cooling time, chosen to be so short that the 
upper boundary stays isothermal, f{z) = (z — zi) / {zq — zi) 
a function which is applied in the interval zq < z < Zi, and 
Co — e(zo) the value of internal energy (temperature) at 
the top of the box. The cooling term serves as a parametri- 
sation of the radiative losses occuring at the surface of the 
star. We assume an ideal gas with 



P = pe(7- 1) 



(10) 



where the ratio of the specific heats 7 — cp/cy — 5/3. eo 
and the stratification in the box are fixed by the parameter 



6 



(7 - l)eo 
gd 



(11) 



which defines the pressure scale height at the top of the 
box. 



2.2. Boundary conditions and initial setup 

We adopt periodic boundary conditions in the horizontal 
directions, and closed stress free boundaries at the top and 
at the bottom 



du,. 



= M, = 



dz dz 
For the vector potential we use 



at z = zq , 23 



(12) 



dA^ _ dAy 



dz 



dz 



= A,=0 



at z = zo,Z3 , (13) 



which constrain the fluctuating magnetic field to be verti- 
cal at the boundaries. We split the heat conduction term 
into two parts 



cond 



= V • [KtV(e - e) + KhVe] , 



(14) 



where the first term acts only on the fluctuations and the 
latter only on the mean, i.e. horizontally averaged, stratifi- 
cation. Thus Kt and can be considered as the turbulent 
and radiative conductivities, which satisfy Kt 2> in real 
stars. We define the conductivies as 



Kt = IPXo , 

(7 - m 



(15) 
(16) 



where xo is the reference value of the thermal diffusivity, 
i^r the input energy fiux, and V the mean logarithmic 
temperature gradient in the initial state. 

As a boundary condition we fix the temperature at the 
top of the box and apply a constant mean temperature 
gradient at the bottom 



eo 



de 



(17) 



where V3 is the logarithmic temperature gradient at the 
lower boundary (see below). 

Initially the radiative flux, J'rad = ^hVe, carries the 
whole energy through the domain. In the present paper 
the uppermost layer is initally isothermal, and the strat- 
ification in the convectively unstable region is described 
by a polytropic index 777.2. Often 7772 = 1 is chosen, but 
this has the disadvantage that the time needed to reach 
the final thermally relaxed state which is much closer to 
the adiabatic stratification becomes large (of the order of 
trfiax = p^fP/nh ~ 4 • lO"^ ^yd/g in the present study, see 



^rclax 

also Brandenburg et al. 2005). Thus in practice, only the 
conductivity Kh strictly corresponds to the case m2 = 1 in 
the initial state, but the thermal stratification corresponds 
to polytropic index 7772 = 1.25 which is already close to the 
thermally saturated state. The initial temperature gradi- 
ent in the lower part of the convectively unstable layer 
and the overshoot layer is calculated from 

V(z) = V3 + i{tanh[4(z„, - z)] + 1}AV , (19) 

where V3 is the imposed gradient, and AV = V2 — V3 
the difference between the temperature gradients in the 
convectively unstable layer and the bottom boundary. Zm 
is calculated on the condition that V = Vad = (7 — l)/7 
at z = Z2 in the initial non-convecting state. The density 
stratification is obtained via the equation of hydrostatic 
equilibrium. 

Furthermore, we decrease the input energy flux by 
a factor of roughly ten in comparison to Papers I and 
II. Decreasing even further would be desirable, given 
the fact that the ratio of the input flux in the model 
to the solar energy flux is still F^/Fq w 10^. However, 
from Eq. I|lt)|) it is seen that F^ cx Kh which leads to 
irciax FjT^. Thus Fi- cannot be decreased by a very 
large amount without a prohibitively large increase in the 
needed computational time. We consider the present setup 
to be a good compromise between real stellar conditions 
and the computational resources required. 

2.3. Dimensionless quantities 

We measure length with respect to the depth of the unsta- 
ble layer, d = Z2 — zi, density in units of the initital value 
at the bottom of the convectively unstable layer, po, and 
acceleration in units of the gravitational acceleration g. 
Futhermore, magnetic permeability is measured in units 
of fiQ and entropy in terms of cp. From these choices it 
follows that time is measured in units of the free fall time, 
yfdfg^ velocity in units of \fdg^ magnetic field in terms of 
\/dpoJIog, and temperature in terms of gd/cp. 

The calculations are controlled by the following dimen- 
sionless parameters. The relative strengths of the thermal 
and magnetic diffusion against kinetic diffusion are mea- 
sured by the kinetic and magnetic Prandtl numbers 

Pr = ^ , (20) 
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where xo is the reference value of the thermal diffusivity 
taken from the middle of the unstably stratified layer. 
Rotation is measured by the Taylor number 

A related dimensionless quantity is the Coriolis number, 
which is the inverse of the Rossby number, Co = 2 fir, 
where r = d/wrms is an estimate of the convective turnover 
time realised in the calculation. 

Convection efficiency is measured by the Rayleigh 
number 

Ra=^, (23) 

where (5 = V — Va is the superadiabaticity, measured as 
the difference between the radiative and the adiabatic log- 
arithmic temperature gradients, and i?ph is the pressure 
scale height, both evaluated from the non-convecting ini- 
tial state in the middle of the unstably stratified layer. 

The strength of the imposed magnetic field is described 
by the Chandrasekhar number 



Ch = 



(24) 



where Bq is the magnitude of the imposed field. 

We measure the strengths of the advection terms ver- 
sus the corresponding diffusion terms in the momentum 
and induction equations by the Reynolds numbers 

Re = , (25) 



Rni 



V 



(26) 



where Urms is the rms value of the velocity fluctuations in 
the convectively unstable layer. The parameters Pr, Pm, 
Ta, Ra, and Ch are used as inputs to the model, whereas 
Re, Rm, and Co are values computed from the actual runs. 

The numerical method used is a modified version 
of that presented in Caunt & Korpi H2001|l . employing 
sixth-order accurate explicit finite differences and a third- 
order accurate Adams-Bashforth-Moulton time stepping 
routine. The code is parallehsed using message passing 
interface (MPI). The calculations were carried out on 
the 34-processor BAGDAD Linux cluster hosted by the 
Kiepenheuer-Institut fiir Sonnenphysik. 

3. The Oij-tensor and methods of analysis 

3.1. The aij -tensor 

In the present study we concentrate only on the first term 
of the expansion (Eq. ^| 

£i — aijBj , (27) 

which, in the present Cartesian geometry, can be written 
in matrix form 



(28) 









a-xz ^ 






(XyX 


ayy 








\ azx 


dzy 


dzz J 





Table 1. Summary of the calculations. The following pa- 
rameters are common for all runs: v = rj = 2 ■ 10^^-\/ gd? ., 
F, = 3_HrVo(.9d)^/', V3 = 0.375, Co = 0.2, = 
lO-'^Vrf/Oom, Pr = 0.4, Pm = 1, (5 w 0.0432, and 
i/ph ~ 0.422d. These choices result in Ch = 0.02 and 
Ra K. 1.03 • 10^. A grid of 128^ is used. The quantities Re, 
Rm, Co and Urms are averages over the convective layer 
and time, whereas Ta and are fixed. We estimate that 
rjTc/ll « 10"'^ in all of the calculations, see Sect. 14.41 



Run 


Re, Rm 


Ta 


Co 




e 


Col-90 


257 


6.5 


10" 


1.02 


0.051 


-90° 


Col-60 


260 


6.5 


10* 


1.00 


0.052 


-60° 


Col-30 


264 


6.5 


10* 


0.99 


0.053 


-30° 


Col-00 


253 


6.5 


10* 


1.02 


0.051 


0° 


Co4-90 


241 


1.0 


10® 


4.40 


0.048 


-90° 


Co4-60 


247 


1.0 


10® 


4.38 


0.049 


-60° 


Co4-30 


267 


1.0 


10® 


4.00 


0.054 


-30° 


Co4-00 


319 


1.0 


10® 


3.30 


0.064 


0° 


ColO-90 


209 


6.5 


10® 


13.2 


0.042 


-90° 


ColO-75 


211 


6.5 


10® 


13.0 


0.042 


-75° 


ColO-60 


213 


6.5 


10® 


12.8 


0.043 


-60° 


ColO-45 


222 


6.5 


10® 


12.3 


0.044 


-45° 


ColO-30 


250 


6.5 


10® 


10.7 


0.050 


-30° 


ColO-15 


283 


6.5 


10® 


9.38 


0.057 


-15° 


ColO-00 


408 


6.5 


10® 


7.23 


0.082 


0° 



where we note that the horizontal mean magnetic fields 
can vary as function of depth due to the contributions 
generated by the fiuid motions in the course of the calcu- 
lation. R suffices to make three calculations with uniform 
magnetic fields imposed along each coordinate axis to ob- 
tain all nine components of the -tensor. This is done 
separately for each rotation rate and latitude. 

Separating the symmetric and antisymmetric parts of 
the tensor one arrives at (see e.g. Radler [l980|l 



(ay + aji)/2 , 
~eijkajk/2 , 



(29) 
(30) 



where aij is a symmetric tensor whose diagonal compo- 
nents describe the generation of the mean magnetic field 
and the off-diagonals contribute to the pumping effect. 
7i encompasses the antisymmetric contribution and de- 
scribes the part of the pumping effect analogous to the 
advection term in the induction equation. More explicitly, 
the electromotive force can now be written as 



£ = aB + 'y X B 



(31) 



where a is a symmetric tensor of second rank and 7 the 
vector defined in Eq. 1)301) . The impact of the off-diagonal 
components of a on the pumping effect can most easily 
be seen by writing (Paper II) 



n 



B 



7 



(x) 



X B, 



- 7*^) X By 



(32) 
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where Bj = BiBj, and 



n 



7 



7' 



7 




(33) 
(34) 
(35) 
(36) 



Eqs. 1(331) to ^ show that the pumping velocity 
can be different for different magnetic field components 
(Kichatinov I1991L Paper II). Note that in contrast to 
Eq. H31|) . which is independent of the chosen coordinate 
system, the coefficients and f^'^^ in the coordinate de- 
pendent Eq. (|32|l are not vectors on account of their trans- 
formation properties. 

In what follows, we refer to the a-effect as the diagonal 
part of aij defined through Eq. H33() , in distinction to many 
earlier studies where it is considered as the full aB con- 
tribution to Eq. (jSU (e.g. Krause & Radler 1980; Radler 
HMH Radler & StepanovEOnSS)- Furthermore, 7, defined 
as the off-diagonal part of a^- via Eq. (|30|l . is from now 
on referred to as the general pumping effect, whereas the 
7*^*-', defined via Eqs. (|34(I - H36|I . describe the field-direction 
dependent pumping effect. 

3.2. Summary of the computations and averaging 
tecliniques 

Table ^ summarises the calculations. A parent run was 
evolved without rotation or magnetic fields until convec- 
tion achieved a statistically steady state. In the present 
case the length of the parent run is i = 100 in units of 
y/d/g. The final state of this run was used as an initial 
condition for the runs in Table ^ for which an appropriate 
rotation vector was imposed. These purely hydrodynam- 
ical runs were evolved up to to — 500-\/ d/g in order to 
minimize the effects of the transients due to the added 
rotation and to allow time for thermal relaxation. After 
this an uniform magnetic field was added. 

Similarly as in Paper II we reset the fiuctuating com- 
ponent of the magnetic field after every T — 50-\/ d/ g time 
units in order to minimize the contamination caused by 
the generated gradients of mean fields. The resetting inter- 
val is two times longer than in Paper II, which is justified 
by the fact that convection, i.e. the average rms velocities, 
in the present calculations are slower by approximately 
the same factor (see Table due to the smaller input 
energy flux F-^. In the averaging, the first At{ = lOy^d/g 
and the last A^i — 30-^/ d/g time units of each subinterval 
are neglected in accordance with the procedure of Paper 



0.; 
o.c 
-0.2 

-O.f 
-0.8 
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Fig. 1. Variation of axx ~ £x/Bx{z) in units of O.Olvd^ 
in the run Col-90 with an imposed field in the x-direction. 
Note that B^(z) = 0, max(|Sj,(z)|/|B^(2)|) = 0(0.1), and 

max{\ ^^^^^^^ \/\Bx{z)\) = O{0.1d-^), indicating that theh 
influence on axx is only minor. The dotted curves show 
the average values from each of the 10 individual subin- 
tervals. The thick solid curve shows the average of the 10 
curves, and the dashed lines indicate the error estimates 
computed from the standard deviation a of the subaver- 
ages according to a/^/N. The thin solid lines indicate the 
full range, denoted from now on by ct, of values attained 
in the iV = 10 subintervals. The vertical dashed lines at 
z — and z = I denote the top and bottom of the con- 
vectively unstable region, respectively. 



II. The generated mean fields in this subinterval are also 
taken into account when computing the coefficients . 
As demonstrated in Paper II, only a few subintervals are 
needed for convergence. In the present study the results 
are averages over 10 subintervals. 

The horizontal and time average of a quantity / over 
the ith subinterval is defined via 



1 



AtL^ 



/ / f{x,y,z,t)dxdydt , (37) 

-iLJ-iLJto+iT+Atf 



where L = Lx — Ly, and At — T ~ Atf — Ati is the length 
of the subinterval in time. Thus we obtain N realisations 
of the quantity /j(z) that can be considered independent. 
An additional average over these N values is simply 



N-l 



/(^) 



(38) 



where f{z) will be the final result. An error estimate can 
be made by computing the standard deviation a of the 
quantity f{z) and dividing it by the square root of the 
number of the independent realisations, or to show the 
full range of values attained in the N subintervals. See 
Fig. n for a typical result. 
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Co = 1 Co = 4 Co = 10 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 



z z z 

Fig. 2. Kinetic helicity H — uj ■ u in units of 0.01 g for the approximate Coriolis numbers 1 (left), 4 (middle), and 10 
(right). Each curve corresponds to a different latitude as indicated in the rightmost panel. The top and bottom of the 
convectively unstable regions are denoted by dashed vertical lines. 



4. Results 

4.1. Kinetic helicity 

Earlier numerical studies of the a-effect have shown that 
the FOSA result, Eq. ||2Jl, which states that a is propor- 
tional to the negative of the kinetic helicity holds at least 
qualitatively (Brandenburg et al. il99G., Paper I). Recent 
local helioseimic studies have been able to extract informa- 
tion on the kinetic helicity, Ji, in the surface layers of the 
Sun using the correlation of the horizontal components of 
the divergence and the curl of the flow as a proxy (Gizon & 
Duvall 1^003(1 . The result is that for the observed latitude 
range up to the latitude Q k, ±45°, the correlation follows 
a cos^ latitude profile. Using numerical convection mod- 
els Egorov et al. H2U04|I showed that the correlation can 
indeed be used to map the helicity. These authors were 
also able to reproduce the observed result with a Coriolis 
number Co w 0.1 which is consistent with the fact that 
the observations correspond to the uppermost 20 Mm of 
the solar surface where the Coriolis number is of this order 
of magnitude. 

However, in a recent study Kapyla et al. 1)2004(1 found 
that the kinetic helicity, averaged over the convectively 
unstable layer, changes sign as the Coriolis number is in- 
creased to the range expected in the deep layers of the 
solar convection zone, Co w 10. This result is in apparent 
contradiction to analytical theory (Riidiger et al. 1999) 
and the observational and numerical results mentioned 
above. This discrepancy, however, turns out to be partly 
due to the volume averaging performed by Kapyla et al. 
H2004|l . although some apparently new features are also 
present. Fig. |21 shows Ti. = uj ■ u for the present calcula- 
tions with Coriolis numbers 1, 4, and 10. It is seen that 
whereas H. is mostly positive at all depths for Co = 1, 
there is a negative region in the deep layers of the convec- 
tion zone and the overshoot layer for more rapid rotation. 
The latitude trend detected by Kapyla et al. (, 2004^) is due 
to the fact that the region of negative helicity expands 



towards the equator whilst the positive peak near the top 
shrinks. Furthermore, the values near the top seem to fol- 
low the more familiar cos trend even for Co k, 10. These 
results suggest that if Eq. ((SJ holds, the latitudinal be- 
haviour of the a-effect should mostly change in the deep 
layers and the overshoot region, but not so much in the up- 
per layers of the convection zone. One must bear in mind, 
however, that Eq. Q holds only for the isotropic case, 
and that anisotropic expressions derived under FOSA give 
significantly different results and are in better agreement 
with the numerical results (see Sect. I4.4ll . 

4.2. The a-effect 

The diagonal components of the tensor describe the 
a-effect which is responsible for the generation of large 
scale magnetic fields in solar and stellar dynamo mod- 
els. In other contexts, we note that dynamos can also 
be supported by the off-diagonal components of ol (e.g. 
Radler & SteDanov l2n06b|l or via the ^Ix J fRadler ll9fi9l 
Roberts and W x J-effects (Rogachevskii 

& Kleeorin 120031 |2004p that involve rotation and shear 
flows in conjunction with large scale gradients of the mean 
fields, respectively. The present results for Uxx, c^yy: and 
azz correspond to age, a^^, and a^r in spherical coor- 
dinates, respectively. Although the present study concen- 
trates on the rapid rotation regime, we have also made 
calculations with smaller Coriolis numbers in order to fa- 
cilitate comparison to earlier results. Fig. |3] shows the a.^ 
components for Co = 1. The horizontal components axx 
and ayy show an expected negative value in the convec- 
tion zone, but basically zero in the overshoot. Latitudinal 
behaviour of the horizontal components is consistent with 
a cos 9 profile (see Fig. (SJ and is in accordance with the 
results for the kinetic helicity (Fig.|2). 

This is of interest for the solar dynamo, where the 
poloidal magnetic field is generated from the toroidal one 
mainly through a^^. To lowest order in rotation, this com- 
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Fig. 3. Diagonal components of the a-tensor in units of 
O.Ol-v/dg for the Col set as functions of depth and in steps 
of A6 = 30° from the south pole (top) to the equator 
(bottom). The dashed vertical hnes denote the top (z = 0) 
and bottom (z = 1) of the convectively unstable region. 
The dotted curves show the full range of values attained 
in the 10 subintervals (see Fig. 



ponent can be represented by a^^ = aig-O oc cos 6, where 
g and Cl are the unit vectors of the gradient of turbulence 
intensity (or equivivalently density) , and angular velocity, 
respectively (Krause & Radler J980) . The sign of the verti- 
cal component, azz, is opposite to those of the horizontal 
components (see also Brandenburg et al. 1990 Paper I) 
and the magnitude of this quantity is comparable to the 
horizontal components in the convectively unstable zone 
and clearly larger, and negative, in the overshoot region. 
The latitude profile is consistent with cos 9 in the convec- 
tion zone and more or less constant in the overshoot layer 
at least up to 6 = —30°. 

The results for Co — 4 are very similar to those ob- 
tained for Co K, 2.4 in Paper II, see Fig. 0] Here the note- 
worthy piece of information is the increased anisotropy of 
the as: whereas for Co = 1 the horizontal components 
were both consistent with a cos 6* latitude profile, for Co 
= 4 only Uxx follows this trend, ctyy, on the other hand, 
is basically constant as function of latitude at least up to 
the latitude —30°. Furthermore, a region of opposite sign 



Fig. 4. Diagonal components of the a-tensor for the Co4 
set, otherwise the same as Fig. |3| 



appears in the overshoot region, consistent with the trend 
seen in the kinetic helicity (Fig. [SJ. The vertical compo- 
nent a^z predominantly shows a mixed sign, but there are 
indications of a region with negative sign near the top, a 
trend that will emerge much more strongly for more rapid 
rotation (see Fig. [SJ. The latitude trend of azz is similar 
to that of ayy. 

The previous numerical studies of the a-effect have not 
gone far beyond Co = 1, although, using mixing length 
models as a guide, it can be estimated that the Coriolis 
number is of the order of 10 or larger in the deep layers 
of the solar convection zone (Kapyla et al. 2()()")1) }. Fig.jS] 
presents results for rapid rotation, i.e. Co = 10. There are 
some distinct differences to the models with slower rota- 
tion. First of all, Uxx changes its sign between latitudes 
45° and 60° in the convection zone. Note that the figure 
shows the full range of values from the ten subintervals, 
and that the actual error is significantly smaller (see Fig. 
O. The value in the overshoot layer, however, still peaks at 
the pole and diminishes monotonically towards the equa- 
tor, although in a rather steeper fashion than cos 6* (see 
Fig- 13 ■ It is also noteworthy that the magnitude of axx in 
the convection zone becomes small in comparison to that 
of a^y. 
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Perhaps the most interesting change that can be ob- 
served in the rapid rotation regime is the strikingly differ- 
ent latitude dependence of Fig- EI); the value in 
the convection zone now peaks at latitude —30°, whereas 
the value in the overshoot is basically constant as function 
of latitude all the way down to O = —15°. Furthermore, 
the sign of the vertical a is now totally reversed in com- 
parison to the case of Co = 1. The maximum values of azz 
in the convection zone and the overshoot layer also occur 
at intermediate latitudes = —30° ... — 45°. 

We note that the trend of the horizontal a-effect as 
a function of rotation at the pole, where a^x = (^yyj 
is consistent with the behaviour found in Paper I. Also 
the sign change of azz is reproduced. Analytical theory 
(Riidiger & Kichatinov I1993|l predicts that the horizon- 
tal a-effect should saturate to a constant value when Co 
is large enough. This seems to occur for ayy in the con- 
vection zone at and near the pole where the values stay 
more or less constant for all Coriolis numbers. The values 
in the overshoot layer and in the convection zone at lower 
latitudes, however, do not seem to show saturation; ac- 
cording to the study of Riidiger & Kichatinov (1993) ayy 
grows roughly 15 per cent when Co is increased from 4 to 
10, whereas in the present results the increase is signifi- 
cantly more. It cannot be ruled out, however, that signs 
of saturation could already be seen if calculations with 
intermediate Coriolis numbers were made or whether sat- 
uration occurs when rotation is more rapid than in the 
present study. The prediction of Riidiger & Kichatinov 
l|1993|) of a vanishing vertical a-effect for rapid rotation is 
not realised in the present results. 

4.3. Magnetic pumping 

The existence of field direction dependent pumping of the 
mean field in magnetoconvection calculations was estab- 
lished in Paper II. We find that the qualitative results 
of the present study are not significantly different from 
those obtained in Paper II. The main results for the gen- 
eral pumping effect, Eq. H30(> . in the rapid rotation case are 
shown in Fig.|Hl The latitudinal pumping effect, j^, is pos- 
titive, i.e. equatorward in the bulk of the convection zone 
and essentially zero in the overshoot layer. The maximum 
of is near the equator, O = —15°, where the value can 
be as high as 25 per cent of the turbulent rms-velocity. 
If we consider the model to describe the deep layers of 
the solar convection zone where u ~ 10ms~^ according 
to mixing length models, the latitudinal pumping velocity 
would correspond to a velocity of 2.5 ms~^, which is of 
the same order of magnitude as the expected values for 
the equatorward meridional flow in the deep layers (e.g. 
Rempel [2005|) . Taking into account the off-diagonal com- 
ponents of the a-tcnsor (Fig. |51) , it is seen that the equa- 

(v) 

forward pumping of the toroidal field, jx — ly ~^ ctyz, 
is further enhanced by the positive ayz- The same term 
enters the equation for the latitudinal pumping of the ver- 

fz) 

tical field with a negative sign and causes 7i to be some- 



what smaller, but due to the small magnitude of ayz in 
comparison to jx, the qualitative results do not change. 
We note that in comparison to Paper II the equatorward 
trend of the general latitudinal pumping effect jx is more 
pronounced in the rapid rotation regime. 

The longitudinal component 7^ is predominantly retro- 
grade in the bulk of the convection zone but often exhibits 
a thin prograde region near the top boundary. The maxi- 
mum occurs again rather near the equator at Q = —15°. 
The latitude and depth distribution of axz is virtually 
identical to jy but with different sign. This causes the lon- 

(x) 

gitudinal pumping of the latitudinal field, 7y , to become 
very small at virtually all latitudes, whereas the vertical 
field is pumped in the retrograde (prograde) direction in 
the convectively unstable (overshoot) layer with a velocity 
comparable to that of latitudinal pumping velocity. 

The vertical pumping effect is directed mainly down- 
wards, as has been found in various earlier numerical 
studies (Nordlund et al. HM^ Brandenburg et al. HM^ 
Tobias et al.^98, 2001J Paper II; Ziegler & RiidigerHISI). 
As was found in Paper II the vertical pumping is domi- 
nated by the diamagnetic effect, which tends to transport 
the magnetic field from regions of large turbulent inten- 
sity, characterised by large rms-value of the fluctuating 
velocity, towards regions of lower intensity (e.g. Radler 
1968; Krause & Radler :198Q'] . In the anisotropic case (see 
Eas. 1471 and I48|) the vertical pumping is proportional to 
—dzul. There is, however, a much stronger latitude de- 
pendence of 7^ than in the more moderate rotation cases 
studied in Paper II, with even a region of upward pump- 
ing in the upper layers of the convection zone at low lat- 
itudes |6| < 15°. The vertical pumping of the latitudinal 
field, jz — Jx + ctxy, is almost constant as function of 
latitude in contrast to the strongly varying general ver- 
tical pumping 7^. The pumping of the azimuthal field, 
7!^'' = "fx ~ oixy, on the other hand, is directed upward at 
latitudes |e| < 45°. 

4.4. Comparison to FOSA 

The exact determination of the transport coefficients, such 
as ciij and bijk, is the most fundamental difficulty in dy- 
namo theory. The problem is basically to choose a closure 
method that is both practical and accurate. This problem 
applies already to Eq. , where usually only the two first 
terms are taken into account. Some indications already 
exist for the need of higher order derivatives in order to 
accurately describe the electromotive force; see e.g. the nu- 
merical geodynamo experiments of Schrinner et al. H2()()5|) . 
The second part of the closure problem is connected with 
the derivation of £ from first principles. The most widely 
used approximation is to consider only the induction equa- 
tion of the small-scale field and neglect all higher than sec- 
ond order fluctuations, which is the first-order smoothing 
approximation (see below). It is possible to extend this 
approach to take into account higher order contributions 
(Nicklaus & Stix ll988|) using the method of a cumulative 
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series expansion (van Kampen ll974alG)l . However, in this 
approach the backreaction of the magnetic field to the flow 
is neglected or at least implicit. Recently, a new approach 
has been introduced (Blackman & Field 2002, see also 
Brandenburg & Subramanian I2005a|l . where instead of 5 
itself, the time derivative £ is computed and the higher 
order contributions, T, are retained via a damping term 
T = —S/t. This approach is often called the 'minimal 
tau-approximation' (hereafter MTA), which leads to ex- 
pressions closely resembling the FOSA results, but with 
the important distinction that due to the fact that the 
equation of motion is also used in the derivation, in ad- 
dition to the kinetic helicity also the small-scale current 
helicity enters the equation for a. Recent numerical stud- 
ies indicate that the current helicity is indeed important 
for the Qf-effect when the magnetic field is dynamically 
important (Brandenburg & Subramanian I2n0,?b|l . 

In order to compare the present numerical results di- 
rectly to the mean-field theory we derive analytical ex- 
pressions describing the coefficients using the first- 
order smoothing approximation (FOSA). Thus we gen- 
eralise Eq. |5J to the case of anisotropic turbulence. The 
reason to compare to FOSA is that it is relatively straigh- 
forward to apply and that the present calculations are in 
the kinematic regime where the influence of the magnetic 
field on the flow is weak and thus the contribution of the 
current helicity is negligible. Furthermore, this compar- 
ison will be useful in evaluating the applicability of the 
approximation itself in dynamo theory. We start from the 
equation of the fluctuating magnetic field 

i} = V X (U X b + u xB + G - TjV X b) , (39) 

where 



G = uxb uxb . (40) 

We drop the term proportional to U for the time being. 
Furthermore, the diffusion term in Eq. (|39|l can be ne- 
glected since "qTc/ll ~ , where we have used Tc — 
2>^fdfg and Ic = lnjlii k, d (see below). G vanishes under 
the assumptions of FOSA, so we are left with 

bi = (SilS-jm - SimS-jl)dj{uiBrn) ■ (41) 

The emf is now calculated from 

£i = £ijkUjbk , (42) 
where 

hk= f hdt. (43) 

Assuming that the correlation time of the turbulence is 
short and that the mean magnetic field B varies slowly, 
Eq. (|43(l can be integrated. Following this procedure the 
ttij tensor components can be derived. Using the assump- 
tion that Uk{t)diUm{t + t')dt' — Uk{t)diUm{t)Tc, with 



Tc having no dependence on fc, Z or m, we find 



axx = o^xx = -Tc{uzdxUy - UydxUz) (44) 

(^yy = %?/ = -Tc{uxdyUz - u^dyUx) (45) 

azz = (yzz = -TciuydzUx - Uxd^Uy) (46) 

axy = Tc{UydyUz + UzdxUx + UzdzUz) (47) 

ayx = -Tc{uzdyUy + UxdxUz + UzdzUz) (48) 

ttxZ = -Tc{UydxUx + UzdzUy) (49) 

azx = Tc{UxdxUy + UydzUz) (50) 

Clyz = Tc{u^dzUx + UxdyUy) (5l) 

Clzy = -Tc{UxdzUz + UydyUx) (52) 



where Tc is the correlation time. In Eqs. (|49|l to (|52|l . terms 
of the form UkdxUk and UkdyUk have vanished due to the 
horizontal periodicity. In practice, we use Tc as a free pa- 
rameter when comparing with the numerical data. 

Figs. EH and ^3 show all nine components of the 
tensor for the calculations Col-30 and ColO-30. We use 
Tc — 3y^d/g in all of the plots and neglect the possible z- 
dependency of Tc. The FOSA expressions for the diagonal 
components of aij seem to best fit to the full numerical 
results of the corresponding quantities, although the fits 
seem to become poorer when rotation is increased. The off- 
diagonal components axy and ayx are reproduced best in 
sign and magnitude, whereas the sign, but not usually the 
magnitude, of axz, o,zx, and azy is in most cases correct. 
For ayz the correspondence is basically non-existent or in 
the few succesful cases it can be considered coincidental. 
The magnitudes of the off-diagonal components already 
show that a universal correlation time that would fit all 
the data cannot be assigned. We note here that the trace 
of aij is proportional to the the negative of the kinetic 
helicity, i.e. 

6ijaij — —T^u} ■ u . (53) 

This relation is compared with numerical results in Fig. 
IT!?] The correspondence is rather good so that at least the 
sign is mostly correct. However, the large peak of the he- 
licity near the surface of the convectively unstable region 
in the moderately and rapidly rotating cases is not well 
reproduced by the trace of aij . This discrepancy, however, 
seems to diminish at lower latitudes (see the rightmost 
panels of Fig. lT^ . 

Although the numerical results differ from the FOSA 
expressions in details, the correspondence is still remark- 
ably good. Very similar results were found also for higher 
Reynolds numbers, with Rm « 400 and 580 from prelim- 
inary higher resolution calculations. It is surprising that 
the same correlation time gives a reasonably good fit for 
the a-effect for all Coriolis numbers, although increasing 
Co is accompanied by a significantly decreased correlation 
time of the convective motions (see Kapyla et al. I2006|l 
and a similar decrease is also seen when MTA is applied 
to the Reynolds stresses (Kapyla et al. I2005a|l . We can 
also estimate the Strouhal number from 

St = UrmsfcfTc , (54) 
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where fcf is the wavenumber where kinetic energy spec- 
trums peaks. The exact value of kf is difficult to assess 
accurately since it tends to vary in time and also in depth. 
However, we note that the power always peaks in rather 
large scales^, i.e. kf — {5 . . . 10) c?~^, with which we obtain 
values of St which are of the order of unity or larger. As 
in Kapyla et al. H2()()f)|l we find that Mrms decreases and 
fcf increases as function of rotation, effectively balancing 
each other so that St is more or less constant as function 
of Co. The values of the Strouhal number obtained by us- 
ing the correlation time from the fits to analytic FOSA 
expressions therefore suggest that higher than second or- 
der terms should be important for the dynamo coefficients 
already in the kinematic regime. 

5. Conclusions 

We use three-dimensional modelling of convection in rect- 
angular boxes at different latitudes and Coriolis numbers 
to study the a-effect and turbulent magnetic pumping. 
The results for the a-effect at moderate and intermediate 
rotation can be summarised as follows 

— As in previous studies (Papers I and II) we find that 
the horizontal components axx and ayy are negative in 
the convection zone on the southern hemisphere and 
follow a more or less cos6' latitude profile for moder- 
ate rotation (Co =1). The latitude trend is in accor- 
dance with the similar distribution of the kinetic helic- 
ity, which, in the simpler case of isotropic turbulence 
is proportional to the negative of the (scalar) a-effect. 
Values of the horizontal as in the overshoot region are 
consistent with zero. 

— For Co = 4 we find that the latitude trend for axx 
remains effectively the same as for Co = 1, whereas 
the magnitude of the quantity increases. For the same 
Coriolis number, the component ayy is more or less 
constant as function of latitude at least up to 8 ~ 
—30°. A pronounced region of oppositely signed hori- 
zontal as appear in the overshoot region with latitude 
distribution similar to the corresponding quantity in 
the convection zone. 

— The vertical component azz is positive in the convec- 
tion zone and negative in the overshoot for Co = 1 and 
exhibits more sign changes as a function of depth for 
Co = 4. 

As noted recently by Kapyla et al. H2004|l . the qualitative 
behaviour of kinetic helicity changes significantly when 
rotation is rapid enough, i.e. Co « 10, raising the question 
of the a-effect in this regime. This regime also coincides 
with the deep layers of the solar convection zone. New 
results for the rapid rotation regime include (see also Fig. 

Ell. 

— The axx component shows diminishing magnitude and 
a sign change in the convection zone near the equator 

^ The smallest possible horizontal wavenumber that fits in 
the box is A; = 2-n jL^ = 7r/2. 



for Co — 10 and a steeper than cos 6* trend in the over- 
shoot. 

— In the stellar context the most interesting component 
is ayy, corresponding to a^^ in spherical coordinates. 
Unlike the earlier studies where rotation is not so 
strong, ayy in the convection zone no longer peaks 
at the poles, but rather much closer to the equator 
at = —30°. The value in the overshoot is virtually 
constant as function of latitude. 

— Another new feature is the sign change of azz , which 
now coincides with the azimuthal component. Also the 
latitudinal distribution is similar, convection zone and 
overshoot values peaking around 8 = —30° ... — 45°. 

— The maximum values of the diagonal components of 

of the order of 0.2 — 0.3Mrms- Considering the 
models to describe the deep layers of the solar convec- 
tion zone where u w 10ms~^, the magnitude of these 
terms would be of the order of 2 — 3 ms~^. 

The turbulent pumping effects can be summarised as fol- 
lows 

— The latitudinal pumping is equatorward and peaks 
near the equator, i.e. around 8 = —15°. The mag- 
nitude of jx is maximally about 0.2Mrms, which would 
corresponds to 2 — 3ms~"'^ in the deep layers of the 
solar convection zone. 

— Longitudinal pumping is retrograde in the convection 
zone and prograde near the upper boundary of the 
domain. At the equator the prograde region is larger. 

— Vertical pumping is predominantly directed down- 
wards and dominated by the diamagnetic effect. 
Whereas the poloidal field is always pumped down- 
wards, the azimuthal field is pumped upwards at low 
latitudes, |8| < 45°. 

The present new results are likely to have consequences 
for mean-field models of the solar dynamo. First of all, 
a long-standing problem of solar dynamo models is that 
the radial differential rotation is strongest at high lati- 
tudes (e.g. Schou et al. 1992|) and the often assumed sim- 
ple form of the a-effect, proportional to cos^, imply that 
the dynamo is most efficient at the poles, but virtually no 
sunspots are observed at latitudes higher than |8| > 40°. 
The tendency of a^y to peak at latitudes around |8| = 30° 
in the regime where the rotational influence is comparable 
to the bottom of the solar convection zone is likely to al- 
leviate this problem. The second major problem for solar 
dynamo models is due to the positive radial gradient of the 
angular velocity near the equator (e.g. Schou et al. I1998|l 
which, with a positive a-effect in the northern hemisphere, 
leads to poleward migration of the activity belts (Parker 
I1987fl . in contradiction with observations. The latitudinal 
pumping will most likely also help in alleviating this prob- 
lem as well, and it has been noted that downward pumping 
has a similar effect in mean-field models (Brandenburg et 
al. I1992(l . The longitudinal pumping velocity 7^ plays no 
role in axissymmetric dynamo models, although its gradi- 
ents and Ve70 give rise to terms analogous to those 
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that appear due to differential rotation. The magnitude 
of this effect, however, is hkely to be small in compari- 
son to that of differential rotation. Longitudinal pumping 
may be important in more general configurations where, 
for example, the propagation of non-axissymmetric struc- 
tures could be explained by this effect. The motivation for 
this is that in the case of the Sun evidence for relatively 
short-lived nonaxisymmetric structures exists (of the or- 
der of 10 rotations from sunspot statistics by Pelt et al. 
H2006|l up to the length of the sunspot cycle of 11 years 
from solar flares by Bai l2003ll . whilst for rapidly rotating 
late-type stars the concept of active longitudes persist- 
ing over several stellar cycles is widely accepted (see e.g. 
Berdyugina & Tuominen 1 1 998|l . 

In order to investigate the applicability of the first- 
order smoothing approximation (FOSA) we have derived 
analytical expressions of the transport coefficients which 
we compare to the quantities obtained numerically. We 
find rather good correspondence for the a-effect in the 
moderate rotation case (Co = 1), whereas the off-diagonal 
components of are only in qualitative accordance or fail 
to reproduce even that. For more rapid rotation the fits 
are poorer, but still manage to capture the qualitative fea- 
tures correctly. The correlation time, used here as a free 
parameter to fit the numerical data to the analytical ex- 
pressions, turns out not to have a universal value. Despite 
of this, we have used Tc = 3^/ d/ g for the whole parameter 
range investigated, as it seems to reproduce the a-effect 
satisfactorily for all rotation rates investigated. With the 
correlation times obtained from the fitting procedure we 
can also estimate the Strouhal number, for which we found 
values of the order of unity independent of rotation. The 
lack of dependence on rotation results in directly from the 
universal correlation time used in the fitting which, how- 
ever, is not consistent with the correlation time computed 
directly from the flow (Kapyla et al. l2005allM7H|l . At the 
moment the reason for this behaviour is unclear. 

In the present study the a-quenching (e.g. Cattaneo & 
Vainshtein I19H1I Vainshtein & Cattaneo I1992|l issue has 
not received any attention, but we note that exploration 
of this issue is further required. A related problem con- 
cerns the magnetic helicity fluxes (Kleeorin et al. 120001 
Vishniac & Cho 20Q1J and the effects of boundary con- 
ditions (Brandenburg & Sandin l^004|l on the efficiency of 
the a-effect. Furthermore, to date, the &ijfc-tensor com- 
ponents have not to our knowledge been calculated from 
numerical models of convection. These points, however, 
cannot be addressed within the scope of the present study 
but should be explored in the future. 
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Fig. 5. Diagonal components of the a-tensor for the ColO 
set in steps of AB — 15° from the south pole (top) to the 
equator (bottom). Otherwise the same as Fig. [3 
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Fig. 6. Diagonal components of the a-tensor for the Col 
set as functions of colatitudc in units of 0. 01-^/35 from the 
top (left panels), and bottom (right panels) of the convec- 
tively unstable region (see Fig. ^ . The error bars denote 
conservative error estimates calculated from the full range 
a, denoted by the dotted lines in Fig. 13 divided by vW, 
where N = 10. The dashed lines show cos 0-latitude pro- 
files for orientation. 
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Fig. 7. Diagonal components of the a-tensor for the ColO 
set as functions of colatitude in units of O.Oly/dg from the 
top (left panels), and bottom (right panels) of the con- 
vectively unstable region (see Fig. [SJ. The dashed curves 
show the expressions denoted in the legends of each panel 
for orientation. Notice the broader plotting range in the 
panels in the lowermost line. The error bars are computed 
as in Fig. El 
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14 



Kapyla et al.: Dynamo coefficients III 



o 0.5 
o 0.0 

° -0.5 
-1,0 



0.0 0.5 1.0 1.; 



P 0.5 
o 0,0 
5 -0.5 
-1.0 

-1.5 



o 0.5 
o 0.0 



-1-0 
-1.5 



3,0 0.5 1-0 1-b 



O.Q 0.5 1.0 



o 0,0 
o -0.5 



0,0 0.5 1.0 1.5 




0.0 O.D 1.0 1.: 




0.0 0.5 1.0 



0.0 0.5 1.0 1.5 



0-0 0-5 1-0 1-5 




0,0 0,5 1-0 1,5 



0.0 0.5 1.0 1.5 




0.0 0.5 1-C l.b 




0.0 0.5 1.0 



].0 0.5 1.0 1.5 



0.0 0-5 1,0 1.5 



0.0 0.5 1.0 



0.0 0.5 1.0 1,5 



D.O 0.5 1.0 1.; 




0.0 0.5 1.0 1.5 



D.O 0.5 1.0 1.5 




0.0 0.5 1,0 1,5 






0.0 0.5 1.0 1.5 



0.0 0-5 1-0 1-5 



D.O 0.5 1.0 1.5 




0,0 0,5 1-0 1,5 



Fig. 10. All nine components of the a^j-tensor for the run 
Col-30 (thick solid curves) and the corresponding first 
order smoothing results (thick dashed curves) according 
to Eqs. (031) to (O- Tc = '^^fdfg in all of the plots. The 
thin dotted and dashed lines give the full range of values 
for the full numerical results and the FOSA equivivalents, 
respectively. 
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Fig. 9. OfF-diagonal components of the a-tensor for the 
ColO set in units of O.Olv^- 
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Fig. 11. All nine components of the a^-tensor for the run 
ColO-30 (thick solid curves) and the corresponding first 
order smoothing results (thick dashed curves) according to 
Eqs. to Tc = iy/d/g in all of the plots. Otherwise 
the same as Fig. ^1 
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Fig. 12. The sum of the diagonal components of Uij (solid 
curves), and —TcOJ ■ u (dashed curves) as functions of 
depth. Both quantities are given in units of Q.Ql\/dq. The 
leftmost column gives data at latitude —90°, the middle 
from —60°, and the rightmost from —30°. The results at 
the equator are not shown. The uppermost panel gives re- 
sults for slow (Co = 1), the middle for moderate (Co = 
4), and the lower for rapid (Co = 10) rotation. 



